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en ■ This paper deals with asymptotic expressions of the Mean Time To Failure (MTTF) 

(__,■ and higher moments for large, recursive, and non-repairable systems in the context 

pLn \ of two-terminal reliability. Our aim is to extend the well-known results of the series 

^^ I and parallel cases. We first consider several exactly solvable configurations of iden- 

Y^ • tical components with exponential failure-time distribution functions to illustrate 

I— I. different (logarithmic or power-law) behaviors as the size of the system, indexed 

, \ by an integer n, increases. The general case is then addressed: it provides a simple 

^ ■ interpretation of the origin of the power-law exponent and an efficient asymptotic 

>0 , expression for the total reliability of large, recursive systems. Finally, we assess 

vQ \ the influence of the non-exponential character of the component reliability on the 

f^ ' n-dependence of the MTTF. 
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1 Introduction 



In non-repairable systems, the Mean Time To Failure (MTTF), i.e., the av- 
erage value ( t ) of the occurrence of failure is a key parameter of the cor- 
responding reliability pl2][3] . Analytical expressions of the MTTF have long 
been well known for simple configurations such as series, parallel, and k-ovX- 
oi-n systems, where each element has a reliability described by the exponential 
distribution p{t) = exp(— At) — the asymptotic dependence of the MTTF for 
a total number n of elements is 1/n (series) and Inn (parallel) — or by more 
complex time distributions [Tlf3] . 
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In this work, we show that these results can be extended to recursive, meshed 
systems of arbitrary size, the latter being indexed by an integer n. Our pa- 
per is organized as follows. In Section [21 we briefly survey definitions and 
simple results that apply for series, parallel and fc-out-of-n systems. General 
expressions for the higher moments /im = (f^) and for the associated cumu- 
lants Km are also given for k-out-oi-n systems. Section [3] is devoted to a few 
exactly solvable configurations [^f5|6|7] of identical components with exponen- 
tial failure-time distribution functions, which give rise to different behaviors 
from those of the series-parallels ones as n goes to infinity. These examples 
allow us to derive the asymptotic expansion of the MTTF when n is large, 
and pave the way to the general case addressed in Section HI which can be 
roughly divided into "series-like" and "parallel-like" configurations. Section [5] 
provides a simple, approximate expression for the corresponding total relia- 
bility and contributes to an understanding of the behavior of large systems, 
which are an important issue [S]. We finally investigate in Section E] how the 
asymptotic n-dependence of the MTTF and higher moments is modified by 
taking non-exponential failure-time distribution functions into account. 



2 Results for simple systems 

2. 1 Definitions 

As explained in many textbooks [T|2|3] ). if R{t) is the system's reliability, the 
MTTF is defined by 

oo oo 

MTTF = {t)= f t {-R'{t) dt) = I dt R{t) ; (1) 



The higher moments /im are 

oo oo 

;,„ = ( r ) = / r {-B!{t) dt) = m I dt r-^ R{t) . (2) 



For instance, the standard deviation a is given by o"^ = (t^) — (t)^. The 
moment generating function f{z) is a formal series defined by 

00 / ±m \ 

/(^) = l + E^^^'"; (3) 

fc=i "^' 



by construction, {f^) = /'•""-'(O) (the mth derivative of f{z) taken at z = 0). 
Since /(O) = 1, the cumulant generating function g{z) is subsequently defined 
by 

/(^)=expff; ^zA=expgiz), (4) 

\k=m / 



where k^ is the cumulant of order m; 1^2 = fJ'2 ~ f^f is the variance, ^3 = 
/is — 3/ii /i2 + 2/if, and so on. Since g{z) = ln( e^* ), k^ = g^^^\0). Using the 
same integration by parts as in eq. ([1]) leads to 

f(z) = (e"*) = l + z f dte''R(t) = l + zR(-z), (5) 





where R is the Laplace transform of the reliability. 

Calculations are often performed by replacing R{t) with an exponential func- 
tion, Rexp{t) = exp(— At). We have then 

00 . 







and ( e^* ) = ( 1 — — j . Consequently, 



</(-P)(^) = -ln(l-^') , (7) 



from which 

(m-l)\ 



^(exp) 



A^ 



When p{t) = exp(— At) is the common reliability of the system's elements, the 
total reliability is 



i?(t)=7^(|>(t))=7^ e^"* , (9) 



where TZ(p) is the reliability polynomial; using the change of variable t 
Inp 

— A~' *^ ^^^ 



oo oo 

MTTF= f dtRit) = - I d{Xt)n 



e-^* 



and 



i/f.(p), (10, 





t™) = ^/-(-lnpr"^7^(p). (11) 

^ 



We make full use of eqs. (11 01) and (jTTi) in the next sections. The moment 
generating function also reads 



oo 







2.2 A few results 

• For n elements in series, TZnip) = P"? so that -R„(t) = exp(— nAt), and 
^rprpp(serics) = J^ [T|3]. Equations( ED and (ED become 

(series) ^ ^^^ ^^3) 

Series) ^ (^-1)! _ (14) 

• If the n elements are in parallel, then -Rn(^) = 1 — (1 — exp(— At))", and 

1 " 1 1 1 / 1 \ 

^rprppCparallel) ^ V- _ ^ (^(^ + l) _ ^(l)) ^lnri+C + — + •■■ (15) 

A^^A AV 2n J 

for n large, where '?/'(2;) = ; is the digamma function and C ~ 

dz 
0.577216 is the Euler gamma constant (see eq. 8.362.1 of [9]). In this case, 

the MTTF diverges as n goes to infinity. 



For k-out-oi-n:G systems, the reliability polynomial is [T][3] 

" n' • 

^m(p) = E 7 '^TriP'i^-Pr~'- (16) 



i=k [ri-zy.z\ 



All the moments could be obtained from eq. (fTTl) . However, the formula 
for the cumulants being eventually extremely simple, we focus on these 
parameters. Starting from eq. (fT2|) . we find 



^ ' AJ pi+^/^ ;r^ (n-z)!?! ^ ^ 



z^ n\ r(z-f)r(n + i-z) 

^A^^, (n-z)!z! r(n + l-f) 



z n\ " r(z 



A r (n + 1 - f ) ,t"fc z! 



(17) 



r(n + l-|) (fc-1)! ' 
the last equality is proven by induction. Consequently, 

'" r(A) r(,, + i-f)' 

whence (see eq. 8.363.8 of [9]) 

-""-^^'Ei; (10) 



^ i=k 



Equation (TTOll generalizes the expressions obtained for the series {k = n 
and parallel (A; = 1) cases, and a similar expression for the variance [10] 

When rra > 1, k^ is bounded by the finite value — C("^)) where ({m 

is the Zeta function. 



Let us now investigate more complex systems. 



3 Simple recursive architectures 



We consider in tliis section configurations tliat are not reducible to series- 
parallel ones. These systems are studied in the two-terminal reliability con- 
text, in which the source and the destination are the large, colored nodes in 
Figs. [TIISIISM and[TDl The size of the system is indexed by an integer n, which 
counts the number of elementary building blocks of the whole structure. All 
the nodes are assumed perfect (they never fail), while the reliability of all the 
edges is described by p{t) = exp(— At) (we may omit in the following the ex- 
plicit reference to time and note this reliability]?, for greater generality). These 
recursive architectures have been solved recently. Our aim is to show that the 
associated MTTFs, as well as moments and cumulants of higher order, can 
be calculated exactly and that their quite distinct asymptotic expansions in 
n give good approximations to the exact results. 



3.1 Ka ladder 




Fig. 1. K4 ladder architecture. 

This architecture, displayed in Fig. [H is constituted by the repetition of perfect 
graphs K4 (the fully-connected graph with four nodes). This configuration 
is exactly solvable even when edges and nodes have distinct reliabilities [B]. 
For perfect nodes, and with all edge reliabilities equal to p, the two-terminal 
reliability has the simple form 



7^n>l(p) = «+(p) Q{p) + «-(p) Cip), 



C±ip) 



p 



2 + Ap- Up^ + 13p'- 4/ ± JA{p) 



(20) 
(21) 

(22) 



., 1+p 2 + 2p+10p^-27p^ + 19p^-Ap^ 

a±{p) = —— ± 1= , 

4 4/4(^ 

A{p)=4: + 32/ - 204 p^ + 452/ - 516/ + 329/ - 112/ + 16|:^3) 



The two eigenvalues are displayed in Fig. [2] as functions of p. As n increases, 
the contribution of (^ prevails over that of (^_. Therefore, when n is large. 



7^4/^a+(/C+(/ 



(24) 




Fig. 2. C-t for the K^ ladder. 

because the contribution of C-{p) vanishes exponentially. The approximation 
given by eq. i^^ lies at the heart of our method for deriving the MTTF's 
asymptotic expansion. Inserting it into eq. flTOj) gives 



t)^\J^a4p)C4p)- 



(25) 



Here, the lower bound (zero) does not play a significant role because C±(0) = 0. 
As n increases, C+{p) is negligible except close top = 1, as illustrated in Fig.O 




Fig. 3. Different powers of C_|_ as a function of p. 



The gist of the asymptotic expansion is therefore to consider a+ and C+ in the 
vicinity of unity. Setting q = 1 — p, we have from eqs. fl2T]) - fl22|) 



a^(l - g) = 1 - 2 g^ + 4 g^ - 3 g^ + 6 g^ + 
-lnC+(l-g) = g^ + 2g5-4g^ + ^gS + --- . 



(26) 

(27) 



Note that C+(l) = cn-(l) = 1 because 7ln{p = 1) = 1. We can write 



At this point, we have to rescale the variable q in order to extract the asymp- 
totic behavior of the integral. Equation ( 1271) gives 

exp[-n (-lnC+(l - q))] = exp[-n(g^ + 2g^ - 4g^ + ■ ■ ■)] ; (29) 



this suggests setting r = nq^, or equivalently q = r^/^n ^/^, so that 

exp[-n(-lnC+(l-g))] = e-^exp[-n(2^^ -4^^ + ■■■)]. (30) 



In the last exponential, the argument is —2 r^'^ n~^''^ + 4 r'^'^ n~^'^ + ■ ■ ■. Be- 
cause of the exp(— r) factor, we can neglect the contribution of large r's, so 
that when n is large, we only need to expand the second exponential and all 
other factors in the limit r — i> 0: 



>>«A/4;^TrTI7l(l-2;^ + 4;;F7I + -)e-'exp|-(2^-4^ + ...)|.(31) 

ni/4 



The error made by replacing the upper bound of the integral by +oo vanishes 
exponentially as n goes to infinity. Keeping only the prevailing terms in each 
factor of eq. (1311) leads to 



1 °° 1 
^ ' X J 4nV4 

= i ^ r(l/4) = i Ml) ^ 0-906402 

For the leading-order term (and this term only), a+ does not play any role 
since it may safely be replaced with 1. The asymptotic n- dependence is not 
1/n or Inn anymore as in the series and parallel cases, but a power-law, n~^^^, 
which slowly decreases with n. The following terms of the expansion may be 
derived easily by expanding all the factors in eq. (1^ : 

r(5/4) 17 r(3/4) _ ^ _ 293 1(5/4) 
^^ nV4 +32 n3/4 4n 512 n^/^ ^"'- ^''''> 



The exact MTTF is obtained straightforwardly by using eq. (fTO!) and the value 
TZn{p) deduced from the three-term recursion relation at the origin of eq. (!20|) : 



T^nip) = iC+ip) + C~iP)) T^n-lip) - C+ip) Up) T^n-2{P) 

= p{2 + Ap- 14p2 + 13p^-4p^)7^„_l(p) 

-p^ (4 - 18p + 36p^ - A2p^ + 30j9^ - 12 p^ + 2/) 7^„_2(p) ,(34) 

with TZoip) = 1 and TZiip) = p{l + 2p — l p^ + 1 p'^ — 2p^). The exact values of 
the MTTF are then obtained by a simple integration of TZn.{p)- The exact and 
the asymptotic (limited to the first three terms of the expansion) results for 
the MTTF are plotted in Fig. HI Even for moderate values of n, the agreement 
between the two is good. 



10 



20 



25 
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Fig. 4. Comparison between exact and asymptotic (see eq. ()33p ) MTTF. 
Following this method, we also find the asymptotic expansion of ( t^ ) by 

and another 1/A in eq. (13T]) . As regards 

rl/4 .1 



adding the factor —2 In 1 



the leading term of the expansion, a mere factor 2 ^^ A~^ is added in the 
integral. Finally 



A2 / 12 s ^ , iL _ 13r(9/4) _ 221 y/^ 
^ ' 2v^ 12 n 10n5/4 480^^/2 



(35) 



Further terms can be routinely obtained using mathematical software such as 
Mathematica. 



After simplification, the variance of the distribution is therefore deduced to 
behave as 



{t')-{t)' 



A2 



TT 21 1 



r(5/4) , ^ , 

2 ^ ' ' ) ^ 12n 



+ iLfl_iZL^U 



16 



.(36) 



We could perform similar calculations for higher moments or cumulants, and 
again would find asymptotic power-law behaviors. 



3.2 Generalized fan 



Sq S^ Sj S3 




Fig. 5. Generalized fan: the source is 5*0, the destination Sn- 

This architecture, displayed in Fig. [5l has been considered in previous studies 
[TT|12fl3j and recently solved for the two-terminal reliability TZn between Sq 
and Sn [5]. For perfect nodes, 



7^. 



p 



p{i-p)y 



+p"(i-p) 



n+2 



np 



1+p' 



(l-p(l-p)) (i-p{i-p)y 



.(37) 



When n ^ 00, TZn clearly tends to the constant TZc 



p 



7^1 



{i-p{i-p)r 

{p (1 — p) is always less than 1/4, so the last contribution in eq. (1371) decreases 
faster than 4~"). This stems from the existence of one path with a finite 
number of hops, namely Sq —* T —* Sn- The MTTF's asymptotic behavior is 
therefore different from that of the preceding section: it does not vary with n. 
This is also true for higher moments, with 



;t-), 



m 
X 



i 

-J dpi-lnpr-' 



p 



{l-p{l-p)) 



2 ' 



(3J 



where {t"^)oo = linin^oo {t"^ )n- The first of these integrals are 



X{t) 



9 + 27rV3 
27 



0.736400 , 



oo = -V^'(l/3) 



4 

27 



TT^^ 0.781302 



(39) 
(40) 



where ip' is the derivative of the digamma function ip. From the first values of 
( t"^ ) , we can infer the general result 

A-(t-)^ = (-l)-^ {1 + ^) (^(-^Hl/3) - ^('"-^ns/S)) 
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ml 



3771—1 



>m-2 



(3 



m.-2 



l)C(m-l). 



(41) 



Depending on the parity of m, the difference ip^^~^\l/3) — ip^"^^^\2/3) may 
actually be further simplified (leaving only ip^'^~^\l/3) for m even, or powers 
of n for m odd). It is easy to prove that, asymptotically. 



'r)c 



ml 



2m Xr. 



(42) 



In that limit, it looks as if only the Sq ^ T ^ Sn connection exists. 



3.3 Double fan 




Fig. 6. Double fan: the source is S, the destination T. 

This configuration, displayed in Fig. [6l is a slight generalization of n double 
links in parallel. As n — > oo, there is an infinity of paths of finite length 
connecting S to T; for perfect nodes, the associated two-terminal reliability 
7^„is[71 



7^. 



l-a+C:-a-C 



(43) 



with 



c 



p 



± ■ 



1,1 

a± = - ± 



1 + 2p (1 - p) ± ^1 + 4p2 (1 _ p)2 



(44) 
(45) 



Here again — if we forget that 1 is a third eigenvalue - we have two eigenvalues 
(±. However, the situation is different from that of the K4 ladder, because 



2 2 ./l + 4p2(i _p)2 



11 



(^_i_ — i> when p ^ 1, while (^ 
that 



1 and a+ ^ 1 when p ^ 0. We also expect 



2\n 



7^n>l-(l-p') 



(46) 



the right-hand side of eq. (146|) corresponding to n elements of reliability p"^ in 
parallel. 




Fig- 7. C± for the double fan. 

As n increases, the contribution of C+ prevails over that of (- (see Fig. [7]), so 
that TZn ~ 1 — «+ C+- Consequently, 



:o.=»i/f(i-..a 



(47) 



Because of the 1/p factor, the asymptotic expansion of ( t )„ is now controlled 
by the behaviors of C+ and «+ ioi p ^ 0: 



C+ = 1 - r - 2p^ + 2p^ + 4p^ - 8/ - ip^ 
a+ = 1 + 2p^ - 8/ + 12/ + 24/ + ■ ■ ■ . 



We can write 



(48) 
(49) 



^2\n 



l-a+C+ = l-(l-r 

+ (1-/)" (l-exp 



-n I 



lnC+ + ln(l -/))]) 



-(l-a+) exp[-n(-lnC+)] 



(50) 



Each term on the right-hand side of eq. (!50|) vanishes for p ^ 0, so that the 
1/p factor does not lead to a diverging integral in eq. (147|) . The first term of 
eq. (1501) gives the prevailing contribution, namely 



12 



1 

A 



dp 
P 



[l-p 



2\n 



1 

2A 

1 
2A 



dr 



ds 



1 " 1 

— Y- 



(51) 



This contribution is — unsurprisingly — half the usual result for n elements in 
parallel, because the reliability p^ translates into a 2 A failure rate. For the two 
other contributions, the change of variable t = np^ gives a factor exp(— r); 
the remaining factors in eq. (150|) must be expanded in the vicinity of r — *> 0, 
as in section 13. 1[ Summing the three contributions gives 



\{t)r 



Inn + C 0r H 95 v^ 1321 
2 ^ 27n ^ 4^^+ i6n3/2 " 24 n^ 



(52) 



A comparison of the exact results with the asymptotic expansion in which we 
have kept the first three terms of eq. fl52|) is plotted in Fig. [HI The agreement 
is good, even for moderate values of n. 




Ill Ti I C 

Fig. 8. Comparison between exact (dots) and asymptotic 1- 

MTTF's for the double fan, in units of A~^. 



vr 
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3.4 Street 3 x n 



C^^^ 9 ^^"^^ ^n 



o p Q g Q T, 



Fig. 9. Street 3 x n: the source is Sq, the destination [/„• 
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In the preceding subsections, we have considered architectures for which TZn 
is exactly known through the analytic expressions of two eigenvalues (±. In 
complex systems, more than two eigenvalues may coexist, and be known only 
as roots of polynomial equations. However, the MTTF's asymptotic expression 
can still be derived from the knowledge of the generating function G{z) [H] 
of the Tori's, namely 

g{z)=Y.n^z^. (53) 



Such is the case of the Street 3 x n, displayed in Fig. [91 This configura- 
tion has been studied for the two-terminal reliability TZn between 5*0 and f/„ 
[6]|15|16|17|18|19|20ll21j (there is actually an offset of 1 between our n and 
these references' n because our source is So). For perfect nodes, Q is given by 
U/iViV2), with [5] 



M = p^ - (1 -p) / (3 + 3p-4p2) z 
+ {l-pfp^ (2 + llp-3p2-2p3) z^ 

+ {l-pfp^ (2-4p + 3p^ + llp^-13j9^ + 3/) z^ 

-(1 - pfp^^ (3 + 6p - 12p2 + lOp^ - 10/ + 4p^) / 

+ {l-pfp^'^ (l + 8p-p2-5p^-p^ + /) z^ 

-(1 - pfp^'' (2 + 5p - 4/) 26 ^ (1 - p)^Vi8 z\ (54) 

Vi = l- (I-p'^Yp [l+p-p^) z 

+ {l-pfp^ {l+p + p^ -2p^) z^ -{l-pfp^z^, (55) 

©2 = 1-^ (2 + 2p + p'^ - Qp^ + bp^) z 

+ (1 -p) p^ (l + 5p + 5p^ -Qp^ -Ibp^ 
+ 13/+/ -2/) z^ 

-{l-pfp'^ (2 + 6p + 6/-26/ + 17/ 



v3 



-18/ + 27/ - 16/ + 3/) z^ 



.4^^5 



+ (1 -pfp^ {l + Qp + Ap^ - / - 17/ 

+9/ + 3/ -2/) z"^ 
-{l-pfp^ (2 + 4p + /-7/ + 3p^ 
+ (1-/'/'^^ (56) 

TV, Vi and V2 are polynomials in both z and p. 
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The eigenvalue of greatest modulus, named C+ again, actually obeys P2 = 
for z = 1/C+ (in the limit p ^ 1, V2 ^ 1 — z and C+ ~^ 1); all other 
eigenvalues tend to zero in that limit. Even though it is not possible to get 
an analytic expression for C_|_ as a function of p {T>2 is of degree 6 in 2;) , we 
can readily compute it numerically. We can also deduce from the constraint 
2^2 (^ = 1/C+) = the expansion of C,+ as a function of q for small g's, starting 
with C+ = 1: 



C+ = 1 - g^ - 4 g^ - 4 g^ + 14 g^ + ■ ■ ■ , (57) 

27 
-lnC+ = g' + 4g^ + 4g5-yg'' + ---. (58) 

The scaling variable r should be n times the leading term of eq. fl58|) . namely 
T = nq^, so that 



= e"^ exp ( - 4 r^/3 n~^/^ - 4 r^/^ ^-2/3 + . . . ^ _ (59) 

a+ is deduced from p and the numerical value of C+ through the residue of Q 
at 2: = 1/C+- The general result is 



where P is the denominator of Q and V = ——. Here, V = 'Di'D2i leading to 

oz 

a+-^ l-2g2-4g=^ + 7g^ + 22g5 + 20g^ + --- . (61) 



After inserting eqs. fl59|) and fl6T|) in eq. fl29l) . and further series expansions 
and integration similar to those of Section 13.11 the final asymptotic expansion 
reads 

, , , r(4/3) 5 r(2/3) 7 , , 

A t„^^^ --^^ + — + ■■•. 62 



Here again, it exhibits a power-law behavior, but n ^/^ this time. In the next 
section, we relate the exponent to a specific property of the network. 
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4 General case 



In the preceding section, we have derived the asymptotic MTTF when the 
two-terminal reliability TZn is known, at least implicitly through a recursion 
relation. Here, we want to show that the leading terms of the MTTF and 
other moments may be obtained for arbitrary, recursive configurations. As 
shown in [B], TZn can be generally expressed as a product of transfer matrices 
— whose size may be large but remains finite. For identical edge reliabilities 
p, the asymptotic behavior is controlled by the largest eigenvalue C+ of the 
(now unique) transfer matrix. We consider in the following architectures that 
look like some kind of "series-like" system, albeit more complex than those of 
the K4 ladder and the Street 3 x n, or to a "parallel-like" one, like the double 
fan. 



4-1 "Series-like" configuration 



This happens when the shortest path connecting the source to the destination 
has a length equivalent to n as n -^ 00. Here, we use again 7^„ ^ 0;+ C+; ^ind 
the MTTF is still controlled by the behavior of C+ and a+ in the vicinity of 
p = 1. The relevant expansions for q ^ have the form 



- In C+(l -q) = aiq' + a^+i g^+^ + ■ 
a+(l - g) = 1 + a'l g + a'2 q^ + 



(63) 
(64) 



from which 



[t), 



1 f dq 



X J 1-q 




'1 + a'iq + a2q^ 



■ ■ ■) exp 



-r2(-lnC+(l -q) - aiq' 



-naiQ' 



.(65) 



The adequate change of variable is now r = naiq\ or equivalently q 



--77; the upper bound of the integral, noj, may again be replaced by 

(nai)V» 

+00. The first term in the asymptotic expansion is therefore 



;On-T 



1 r 1 



X J i (naj)^/' 

1 r(i + i/z) 

A (nai)Vi 



r^/*"irfre"" 



(66) 
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Likewise, in order to calculate the first term of the asymptotic expansion of 
A^ (t^ )„, we have an additional factor —2 ln(l — g), which is equivalent to 2g 
when g — > and brings an extra r^/*. Finally, 



:t')n^i. 



1 r(i + 2/z) 



A2 {naiy/' 



(67) 



so that the variance goes as 



'.t')n-{t)i^^ 



1 r(i + 2/z)-r(i + i/z)^ 



(raaj)^/* 



(68) 



from which 



{t)n 



r(i + 2/z) 
^r(i + i/z)2 



(69) 



which is independent of n. Actually, it only depends on i, which is the lowest 
order of the g-dependence of 1 — C+ when q ^ (see eq. flB5]) ). 

For the higher moments {t'^)n, the generalization is straightforward; we can 
also go beyond the first order in the expansion, following the recipe of the 
preceding section. Setting r] = {nai)~^^\ we find 



A™(n„ = r7^ 



r(i + - 



m 



+r] 



9 m 
+r] — 



"^■^^rfi + i 



OLi 



m 



10 + llm + 3m2 + 12(l + m)a; + 24a^ /2 + m 
24 V^ 

fai+2 ^ fti+i 1 + m + 2 a[\ y [-^ , 2 + m 



2 V a,- 



12 



2 

2 + m 



(70) 



We may wonder: is this result merely formal, or is it actually possible to 
determine A™ (t"*)n for an arbitrary, recursive network ? The answer to this 
question is yes. Even though we do not know the exact value of TZn or the 
associated greatest eigenvalue C+, we can still infer the a^'s and a^s appearing 
in eqs. (I63ll - (1M1) because for large n, TZn ~ ct+Cl- These parameters can be 
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deduced from the expansion of the unavailabihty W„ = 1 — IZn for g — > 0. We 
have 



= 1 - {I + a iq + a' 2q^ + ■ ■ ■) exp \-n (at q' + a^+i g*+^ + ■ ■ ■)! ( ''!) 



and must keep track of the successive powers of g, along with their dependence 
with n. Let us illustrate this claim with the Street 3 x n case. A simple cut 
enumeration gives (when n is large, so as to avoid "boundary" effects) 

Un = 2q^ + in + A) q"" + {An - 7) q^ + ■ ■ ■ . (72) 



The first term of the right-hand side of eq. f[72|) is easy to obtain. There are only 
two cuts of order 2 preventing a connection between source and destination 
(see the top of Fig. [TU]) . hence the 2 q^ term. For the cuts of order 3, different 
possibilities occur as displayed at the bottom of Fig. [TUl Firstly, three parallel 
links may fail; there are n such instances. Secondly, close to the source or the 
destination, there are four triple failures (only two are represented by green 
stars in Fig. [TOl the remaining ones can be deduced by symmetry). This gives 
the {n + 4) q^ term. A comparison between eqs. (17T]) and (1721) gives a'l = 0, 
a' 2 = —2, a' 2, = —4, and 0:3 = 1. From these values, we get 

-<n-;^(r(i.f)-^J^r(i.^)....), (.3) 



which agrees with eq. fl62l) when m = 1. 

In conclusion, the exponent of the power-law behavior in n is nothing but the 
inverse of the number of necessary cuts to isolate each elementary cell from 
its neighbors. Note that Oj is not necessarily equal to 1, since it represents the 
number of independent cuts of order i. 



4-2 "Parallel-like'' configuration 



We can also consider configurations where the reliability is asymptotically 
equal to 1 when n — > 00, as in section 13. 3[ In the vicinity oi p —* (the 
relevant domain here) , 7^„ ~ 1 — a+ C+ and the needed expansions are 



lnC+ = P^p' + P^+lP'^' + ■■■ , (74) 

a+ = l + (3\p + (3'^p' + ---, (75) 






<:^-^<:^^<^^^:^^^ 



o c^ o 










o^ 



o-v-o c^-^<^-^o 




Fig. 10. Cuts of order 2 (top) and 3 (bottom) for the Street 3 x n. Failed edges are 
indicated by stars and crosses. 



so that 



1 — a 



+c: = (i - (i-pT'^) + ((1 -pT'^ - c) + (1 - «+)c: • (76) 



The contribution of A is (omitting the A ^ factor) 



A-/^(l-(l-pr'")=-E^--flnAn + C + -i 
J p ^ ^ I t^^k I \ 2(3, 



n 



(77) 



The contribution of B depends on the value of i because 

1 



lnC+ + Aln(l-p*)=A+iP 



»+i 



A(-2P'^ + 



(75 



If 2 = 1, we must take the two terms of degree 2 into account; otherwise, only 
the /3i+ip*^^ term needs be kept. After asymptotic expansions similar to those 
performed in the preceding sections, we have 



B^ /?2- 



A 



A 



j+i 



r 1 



1 1 

1 



+ ■ 
1 



A 



l+l/i 



n-1/^ + 



(^ = 1) 
f^>l) 



(79) 
(80) 
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The contribution of C is easier to compute because 1 — a+ vanishes as p — *> 0, 
thereby compensating the singular term 1/p in the integral. With the change 
of variable t = n Pip\ we get 

C^ , ~fl.. T(l + -]+---. (81) 

(nA)iA \ ij ^ ^ 



Note that the n-dependence of C is n ^/\ which decreases less rapidly than 
1/n if z > 1. The sum of contributions A, B, and C finally expands as (for 
i> 1) 

A MTTF,. -. 1 (l„(ft „) + C) + (fii - /?;) '-^ + . . . (82) 



As in the preceding subsection 14. 11 the coefficients Pj, /3^, etc. may be deduced 
by evaluating the availability through a path enumeration in the limit p ^ 0. 
For instance, in the case of the double fan 

n^ = np'^ + 2{n-l)p^ + ■■■ (83) 

This must be compatible with the expansions of a+ and — In C+ in eqs. ( TMl) - 
( 1751) . Because eq. (l83l) has no linear term, /5i = /3( = 0. The coefficient of p^ 
being equal to n, we deduce z = 2, /32 = 1, and P2 = 0; the coefficient of p^ 
then implies jSs = jS'^ = 2. Inserting these values in eq. fl82l) gives back the first 
two terms of eq. fl52|) . 



5 Approximate reliability of large, recursive, "series-like" systems 



We have seen that the asymptotic expansion of the MTTF and the higher 
moments for a large, recursive system can be obtained with minimal effort. In 
the same line of thought, is it possible to find an approximate reliability such 
that all its moments give the same value as the true reliability, at least for the 
first terms of the expansion in n. The first term of eq. fITOl) . namely 



reminds us of what would be obtained for a WeibuU distribution [T|3] . Indeed, 
it is straightforward to show that 

i?i°) it) = exp (-n ai X' f) , (85) 
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would give eq. (18^ exactly. 

Is it possible to improve this expression, i.e., propose an effective reliability 
leading to the correct first two terms in the asymptotic expansion of each 
moment {t"^)n ? Provided that a'l = 0, i.e., that there is no cut of order one 
(it would then be easy to factor out this link contribution, and proceed with 
the remaining parts of the system), the answer is again positive. 

Calculating the moments of 



R^:\t) 



exp 



i+l J.^+l^ 



-n [ai A* f + ai+i X'^^ t 



(86) 



we find that 



^"' <*"•>- -(^TT^rO + 7 



a,- 



j+i 



ai [nai]'^/^ \ i 



m f 771 + 1 

— r 1 + ^^- 



.(87) 



Equations (1701) and (1871) match if 
5j+i Oij+i 1 m + 1 + 2 a'l 



ai ai 2 



m+l 



SO that for a'l = 0, the constraint is satisfied when 



fti+i = Oj+i - -ai 



(89) 



This finally gives 



R^\t) 



exp 



-n [a,\'f+[ai+^--ai] A'+4*+i 



(90) 



This expression slightly improves over the Weibull distribution of eq. (1851) . In 
the case of the Street 3 x n, 



R^n\t) 



exp 






(91) 



This expression and the true reliability are plotted in Fig. [TTl the agreement 
is already satisfying for n = 30. 
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Fig. 11. Exact (dots) and asymptotic (full line) reliability for the Street 3 x n 
architecture (n = 30). 

6 Non-exponential distribution functions 



In the preceding sections, we have considered elements whose reliability is 
p{t) = exp(— At). Although this distribution is often chosen because calcula- 
tions are simpler, other models may be used: Weibull, gamma, lognormal, etc. 
pH3] . We investigate here the influence of the true p{t) on the n-dependence 
of the MTTF and higher moments for "series-like" systems. 

We can invert p(t) as t = x(p)) so that eq. (fTOj) transforms into 

1 
MTTF„ ={t)n = J dpTZ^ip) i-x'ip)) ■ (92) 



For higher moments, we would get 



{t^)^ = m j dpx'^'-'ip) n^{p) i-x'ip)) . 





(93) 



For "series-like" configurations, we have again to consider what happens for 
p ^ 1, or equivalently for t ^ 0. Assuming that asymptotically 



-X'(p)^a/3(1-P)^. 



(94) 



we have (because x(l) = 0) 



xip) = t 



(3 + 1 



'1-pY 



(95) 
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so that 



/(/5,l)Ai/(/3+i) 



Keeping — In ^+ -^ ai g* and a+ — > 1 , we get to lowest order 

1 / \ m—l 

( r ). - m / rfg f ^ j g(^+^) (™-^) a, g^ exp {-n a. q^ 

The n-dependence of ( t"^ )„ is therefore affected by the structure of the graph 

(through i and a^) and by the true failure-time distribution of each element 

(through a/3 and /5). The exponent of the power-law asymptotic behavior is 

(P + l)m , ,. 1 .,. .1, 
: , and the total reliability goes asymptotically as 



Rn{t) ^ exp 



-n ai 

V «/3 / 



(9J 



7 Conclusion 



We have shown that very simple asymptotic expansions may be obtained for 
the mean time to failure (and higher moments) for general recursive, meshed 
networks. By contrast with the simple series and parallel systems considered 
in many textbooks, the size-dependence of the MTTF of a "series-like" system 
follows a power-law behavior, whose exponent is linked to the number of cuts 
disconnecting an elementary cell to its neighbors. Comparison with the exact 
results for various architectures show that the agreement is often reached 
when the system contains a few dozens of the (repeated) pattern structure. A 
simple, approximate expression for the effective global reliability of the system 
has also been proposed, which is very simple to derive by a mere enumeration 
of cut-sets or path-sets. 

The calculations have been performed in the context of the two-terminal reli- 
ability of general systems; they obviously apply to all-terminal reliability, and 
would belong to the "series-like" category. 
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